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Long -distance, high-capacity transmission via the TE \ mode in circular 
waveguide is an attractive goal because the theoretical attenuation caused by 
heat loss decreases monotonically as the frequency of operation increases. 
However, continuous random mode conversion caused by the manufacturing 
or the laying of the waveguide causes severe random fluctuations with fre- 
quency in the transfer function of present-day waveguide, a fact which 
apparently eliminates all but the "toughest" modulation schemes, such as 
pulse code modulation and angle modulation. 

In this paper, we present the derivation of a technique for analyzing the 
effects of continuous random mode conversion on an angle -modulated wave. 
The random coupling coefficient is assumed to have a Gaussian probability 
density. The modulating signal is assumed to be a single-sideband frequency 
division midtiplex which may be simidated by a band of random Gaussian 
noise. The quantity of interest is the expected value of the interchannel 
interference noise which appears at baseband. A self-contained section on 
principal residts and examples is presented. The examples use typical 
state-of-the-art data and are significant in that they demonstrate that FM 
appears to be an attractive modulation scheme. They are, however, not 
meant to be an exhaustive study of the problem. Copper waveguide (with 
and without mode filters) and helix waveguide are considered. 

The Fourier series, multiple-echo representation of a stochastic transfer 
function presented in this paper should prove usefid in handling other 
modulation schemes as well as other random channels. 

I. INTRODUCTION 

Long-distance, high-capacity transmission via the TE i mode in the 
millimeter-wave circular-electric waveguide is an attractive goal be- 
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cause the theoretical attenuation caused by heat loss decreases mono- 
tonically as the frequency of operation increases. However, in order to 
obtain sufficiently small heat loss attenuation and delay distortion, the 
operating frequency must be well above the TE i cutoff frequency and 
far into the multimode region. For example, a two-inch I.D. perfect 
copper circular guide will have a theoretical TEoi heat loss of 1.54 db 
per mile at 55 kmc and will propagate 223 additional (spurious) modes. 1 
Real waveguides are never perfect right-circular cylinders but have 
continuous random imperfections arising from the manufacturing or the 
laying of the waveguide. These random imperfections cause conversion 
of the TEoi signal mode to spurious modes. 2 We call the attenuation 
produced by this energy loss "conversion attenuation." Usually more 
important, the imperfections also cause converted energy in the spurious 
modes to be reconverted into the TEoi signal mode. Since the group 
velocities of the various modes will be different, the reconverted energy 
from each mode will be delayed or advanced in time with respect to the 
energy which has remained in the TEoi signal mode. Thus we have a 
signal plus a continuum of echoes. The effect causes the amplitude and 
phase of the TEoi wave to fluctuate randomly with frequency. A transfer 
function which is a random variable can be defined to take this fluctua- 
tion into account in a statistical sense. We call it the "reconversion 
transfer function." 

The fluctuations in the reconversion transfer function distort any 
signal which might be transmitted through the waveguide. A "tough" 
modulation scheme is required in order to overcome this effect. In this 
paper we consider the use of large-index frequency modulation (FM). 
The recent work of Rowe and Waiters 2 provides us with a sound and 
tractable mathematical analysis of continuous random mode conversion. 
We shall use the results of their work as our starting point. The base- 
band signal used to frequency or phase modulate the carrier wave is 
assumed to consist of a single-sideband frequency division multiplex 
of telephone channels. For analysis purposes we can approximate this 
baseband signal by a band of random Gaussian noise. 3 When such an 
angle-modulated wave is passed through a network whose attenuation 
and phase are nonlinear functions of frequency, interchannel interference 
appears in the demodulated signal. In way of explanation, imagine that 
the energy in the baseband modulating signal corresponding to some 
particular telephone channel is removed. When such a signal is impressed 
on a frequency modulator and the resulting wave is transmitted through 
the network in question, detected, and finally demodulated, the received 
output in the originally clear channel represents interchannel interfer- 
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ence. Measurements of this type have been discussed by Albersheim 
and Schafer. 4 

To the author's knowledge, there is only one tractable technique for 
determining interchannel interference in angle-modulated systems which 
does not require that the fluctuation rate of the attenuation and phase 
functions be slow over the bandwidth occupied by the signal. This 
exception is the method of equivalent echoes. 56 Since the phase and 
attenuation of the reconversion transfer function fluctuates rapidly 
across the bandwidths of interest, we are necessarily confined to this 
technique or variations thereof. We see that some sort of multiple-echo 
approximation to the reconversion transfer function is required if any 
accuracy in our results is to be expected. How many echoes are required, 
and how should they be distributed and weighted? The answer is clear 
when we recall that a small echo produces or is equivalent to a sinusoidal 
variation with frequency in the transfer function, i.e., the Fourier trans- 
form of the echo. We simply expand the reconversion transfer function 
in a Fourier series over a frequency interval which we know to be suffi- 
ciently large to contain all of the significant energy in the FM wave. The 
Fourier transform of this series is the desired train of echoes and is an 
exact representation over the bandwidth under consideration. The echoes 
are regularly spaced, and their amplitudes are proportional to the magni- 
tude of the coefficients of the Fourier series. The problem is now reduced 
to obtaining the Fourier series coefficients, and until recently this would 
have been a major stumbling block. Fortunately, however, Howe, 
Warters, and Young have recently made important contributions 2 ' 78 to 
the understanding of continuous random mode conversion. They have 
demonstrated that over the moderate percentage bandwidths in which 
we are interested, the effects of random mode conversion give rise to a 
transfer function which is a random function of frequency much as a 
noise wave is a random function of time. Using their results, we can 
evaluate the statistics of the Fourier series coefficients, i.e., echo ampli- 
tudes, which are, of course, random variables. Moreover, if the cross- 
correlation between coefficients is negligible, we can determine the ex- 
pected value of the interchannel interference for each echo individually 
and then add the results to obtain the total expected value of the inter- 
channel interference. The random coupling coefficient is assumed to 
have a Gaussian probability density. The effects of delay distortion 
caused by waveguide cutoff dispersion are not included in this paper. 
The system is assumed to be equalized. 

Section II presents the principal results of this paper and contains 
examples demonstrating their use. The examples use typical state-of- 
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the-art data obtained from S. E. Miller 9 and are significant in that they 
demonstrate that FM appears to be an attractive modulation scheme 
for use on the millimeter-wave circular-electric waveguide as far as 
continuous random mode conversion is concerned. Both helix waveguide 
and copper waveguide with mode filters are considered. The helix wave- 
guide has superior performance, of course, but copper waveguide with 
periodically placed mode filters may be useful in some situations. 
Section III presents the mathematical development of the results. 

II. PRINCIPAL RESULTS 

The present section contains a statement and discussion of the princi- 
pal results of this paper. We will restrict the discussion to the case of a 
single spurious mode and a single polarization. If other polarizations are 
present, they can be considered individually and the results added on a 
power basis.* The same is true for our purposes for different spurious 
modes. It is certainly true when the spurious modes arise because of 
different and uncorrelated physical imperfections in the waveguide. For 
instance, the TE02 spurious mode is caused primarily by random diam- 
eter variations. Its reconversion transfer function would be uncorre- 
lated with that of the TEu spurious mode, which is caused primarily 
by random straightness variations, if the diameter variations are un- 
correlated with the straightness variations. Even when the spurious 
modes arise because of the same physical imperfections, it can be shown f 
that the cross-correlation in the same frequency interval is negligible. 

There are two cases of interest. In the first, the difference between the 
heat loss attenuation constant of the TE01 and the spurious mode, i.e., 
the differential attenuation constant AA, is uniform. An example would 
be the TE 0n family on both the helix waveguide and copper waveguide 
with helix mode filters. The second case of interest is when the differen- 
tial attenuation is effectively zero, except at periodically located points 
along the waveguide where a total differential attenuation of k is in- 
troduced. An example of this type of behavior would be the TE in family 
on a waveguide having periodically placed helix mode filters. 

2.1 Uniform Differential Attenuation Constant 

Consider a length L of waveguide having a differential attenuation 
constant of AA for the spurious mode of interest. Assume that the ran- 
dom coupling coefficient between the TE i and spurious mode is Gauss- 

* See Ref. 2, Section 4.3.2. 
t See Ref. 2, Appendix G. 
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ian, pure real or pure imaginary, and has a white power spectrum of 
density S e . If the differential attenuation of the total length of wave- 
guide is much larger than one, i.e^, | AA \ L » 1, then the mean-square 
value of the equivalent echoes r» a , normalized by the signal power, is 
given by the expression 

r7= (T/r )(P x /P )r MTho) (1) 

if P x /P «l, where 

P*/Po = reconversion echo-to-signal power ratio (2) 

= S c 2 L/2 | AA | (2) 

T-o - \vl - vq I /4c I Ai4 I (3) 
and 

n = positive integer designating the nth echo. (The echoes lag the 
signal if v 2 > v 2 and lead if v Q 2 > v 2 . Because the interchannel 
interference is independent of this, we use the absolute value | v 2 — v 2 \ 
in (3) and require n to be positive.) 

P x = reconversion echo power, i.e., the power in that portion of the 
signal which has been converted and reconverted once and only once. 

P = signal power which has suffered no conversion. This will hence- 
forth be considered the signal power. 

v , v x = mode cutoff factors* of the TE 01 and spurious modes at 
the band center frequency. This give rise to the differential phase shift. 

c = velocity of light in free space. 

T = 1/Af = time delay between successive echoes. 

Equation (2) has also been derived independently by Miller, 9 using a 
different approach. Equation (1) was developed from a Fourier series 
expansion of the reconversion transfer function over the frequency 
interval Af. Hence, A/ must satisfy certain conditions if meaningful 
answers are to be obtained. First, all of the significant energy in the 
angle-modulated signal must be contained within the bandwidth A/. 
Thus we require a bandwidth of 

A/ ^ 2(f b + 4<r) (4) 

where f b is the baseband bandwidth and a is the root-mean-square fre- 
quency deviation of the angle-modulated wave. This is simply a state- 
ment of the familiar Carson bandwidth. Second, the mathematical 

* See Ref . 2, Appendix A. 
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derivation of ( 1 ) requires the differential phase constant to be approxi- 
mately proportional to frequency over the bandwidth used in the Fourier 
series expansion. This sets an upper limit on A/. In order to satisfy this 
requirement we require that the percentage bandwidth and the square 
of the mode cutoff factors of the modes involved be small compared to 
unity, i.e. 

A//2/o « 1, v<? « 1 and v? « 1. (5) 

This will usually be true for the bandwidth and modes of interest, but 
it should always be checked as a safeguard. Third, the Fourier series 
coefficients must be uncorrected with one another. It is well known that 
Fourier series coefficients become uncorrelated as the expansion interval 
approaches infinity. In a practical case, "infinity" is a number which 
is much larger than the correlation interval of the random variable in 
question. In the present case the autocorrelation function of the recon- 
version transfer function varies as 1/[1 + (2irT f) 2 ], where /' is meas- 
ured about the band center frequency. Thus in order for the cross-cor- 
relation between coefficients to be negligible we require that 



(ttt A/) 2 »1. 



(6) 



Lastly, the expression ( 1 ) for r„ 2 is limited by the accuracy of the per- 
turbation theory from which it was developed. It should be good, how- 
ever, as long as P /P x » 1 as calculated from (2). A plot of the equiva- 
lent echo train is shown in Fig. 1. We see that r is the "time constant" 
associated with the envelope of the echo train, i.e., it is the value of 
time delay where the envelope of the echo train caused by reconversion 
has dropped to the relative value 1/e. 

If the cross-correlation between Fourier series coefficients, and hence 
echoes, is negligible, we may substitute the expression (1) directly 
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Fig. 1 — A plot of the reconversion echo train for uniform differential attenua- 
tion. 
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into the formulas given by Bennett et al. 5 and sum the results to deter- 
mine the total baseband interchannel interference in the top channel. 

As an illustration, for a two-inch I.D. helix waveguide at 55 kmc, the 
dominant spurious mode is the TE 2 , as far as reconversion echoes are 
concerned. Typical data are 9 

| AA | =4.3 db/mile = 0.495 neper/mile 
S c = 0.01 db/mile = 0.00115 neper/mile 
vo 2 = 0.0172 « 1 
v? = 0.058 « 1 
T = 105.8 X 10 -7 sec. 

For a waveguide length of 4000 miles, we have from (2) 

P /P x = 22.6 db»l. 

Let us assume that a signal is applied to this waveguide which has been 
frequency modulated by a rectangular band of Gaussian noise. Let the 
baseband width be f b = 10 mc, and the rms frequency deviation be 
a = 2f b = 20 mc. The significant energy in this wave would be contained 
in a bandwidth of A/ = 10a - 20/ 6 = 200 mc. Thus, if we pick A/ = 
200 mc we also satisfy (5) and (6), i.e. 

A//2/ = 0.002 « 1, (tttqA/) 2 = 4,450 » 1. 

From Fig. 5.7 of Bennett et al. 6 we can read off the factor s(nf b T), 
which multiplies the mean-square value of each echo to give the inter- 
channel interference-to-signal power ratio for that particular echo. 
s(nf b T) takes into account the noise suppressing property of FM as 
well as the variation caused by the magnitude of the delay. A plot of 
s(nf b T) is shown in Fig. 2 for a/f b = 2 and 3. The total interchannel 
interference-to-signal power ratio is given by the sum 

I/S = £ s(nf b T)7? (7) 

if P x /P « 1. We see that since s(nf b T) is essentially constant at —15 
db for nf b T > 1/10, i.e., n > 1, we have to a good approximation 

I/S = s(nf b T) I) ^ 

(8) 

- s(nf b T)(P x /P ). 
In decibels 

I/S = -15 - 22.6 = -37.6 db. 
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Fig. 2 — Weighting function s(nfbT) which multiplies individual echoes to 
give interchannel interference in top baseband channel for FM. 

This corresponds to 27.4 dba at zero level. Thus, other things neglected, 
continuous random mode conversion does not preclude a 4000-mile 
two-inch I.D. helix system. 

If this interchannel interference noise-to-signal ratio were too large, 
one could decrease it still further simply by increasing the modulation 
index. In other words, for the conditions of this example we are able to 
obtain FM improvement against reconversion echoes as well as additive 
thermal noise. However, from Fig. 2 we see that for nf b T « 1/10 we 
would not be able to obtain this improvement. In fact, increasing the 
modulation index degrades the situation in that case. 



2.2 Mode Filters 

Mode filters may be placed at periodic intervals along the waveguide 
to provide large attenuation to the spurious modes. We assume here 
that each filter contributes k nepers of differential attenuation and that 
the waveguide between niters contributes essentially none. Let the 
spacing between mode niters be L and the length of each mode filter 
be negligible in comparison. In the previous section, for uniform differen- 
tial attenuation we found that the envelope of the equivalent echoes, 
when plotted versus delay, was an exponential. In the present case we 
find that the envelope consists of a piecewise approximation to an 
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exponential using straight line segments (see Fig. 3). The exponential 
is e~ * and the straight line segments connect the points defined by 
letting i = 0, I, 2, ■ ■ • , etc. Let us number the trapezoids represented 
by the area under each straight line segment according to the value of i 
at the beginning of each trapezoid. The equation for the straight line 
in the it\\ trapezoid, i.e., between f 2 and € -2(,+I)fc on the ordinate, is 



?/.(r/r 1 ) = «-»*[] + [1 - 6-*][t - (r/n))} 



(9) 



where t x is the length of the base of each trapezoid. The expression for 
the individual echoes is 



2 1 Sc'LLo — 2 ft- 

r„i = e 



[1 +[l -€-*][z- (nr/n)]] 



(10) 



where i is chosen from the inequality i ^ (nT/n) ^ i + 1. Here 
L is the length of waveguide between mode filters, L is the total length 
of waveguide and rj = \ v x 2 — v 2 \ L /2c. All other quantities are as 
defined in the previous section. As in the previous section, we require 

that 



Bb 1 «l, 



v* « 1 , 



A//2/o « 1 and A/ ^ 2(/„ + 4a ) . (11) 



The requirement on &f to ensure negligible cross-correlation between 
the Fourier coefficients of the reconversion transfer function, i.e., be- 
tween individual echoes, is changed, however. The autocorrelation 
function is given by (42). The series converges quite rapidly for large 



-jr /(TS 2 C LL ) 




L,r/r, 



Fig. 3 — For periodically placed mode filters of k nepers differential attenua- 
tion each, the envelope of the equivalent echoes consists of a straight-line ap- 
proximation ?/,• (t/ti) to the exponential «~ 2,fc . 
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differential attenuation per mode filter k. In most cases of interest only 
the first term or two is required. In order for A/ to be much larger than 
the correlation interval we require 

(7TT3A//2) 2 » 1 (12) 

since | f | L = fn . 

Equation (10) requires the number of mode filters N to be large, 
i.e., JV » 1, and that the square of the total differential loss of the N 
filters be large, i.e., e~ 2kN « 1. Under these conditions, the reconversion 
echo-to-signal power ratio is given by 

P x /P = (LL Q S c 2 /2) coth k (13) 

which must be small compared to unity in order for (10) to be accurate. 
As a check, notice that if k « 1, we have 

P 8 /P = LS 2 /2\AA | (14) 

where | AA | = k/L . This is precisely the expression presented in the 
previous section for a constant differential attenuation. 

It is frequently necessary to determine the ratio of the power in the 
ith. trapezoid to the total. From simple geometric considerations this is 
given by 

P X{ /P X = e~ 2ik [1 - r 2fc l (15) 

Equations (13) and (15) have also been derived independently by S. E. 
Miller 9 using a different approach. 

As an illustration, consider a two-inch I.D. copper waveguide system 
at 55 kmc having 15-db mode filters every 200 feet. The TE U , TEi 2 , 
and TMn spurious modes would all contribute significantly because of 
random straightness deviations. The energy in all of these modes will 
be lumped together and considered as if it were in the TE J2 mode. This 
mode has the largest delay and consequently would produce the most 
interchannel interference. Thus we are assured that our calculations 
will be pessimistic. Numerical values were obtained from Miller 

k = 15 db = 1.73 nepers 
U = 200 feet = 1/26.4 miles 
S c = 2 db/mile = 0.23 neper/mile 
vo = 0.0172 « 1 
v 2 = 0.0334 « 1 
T-i = 16.08 X lO" 10 sec. 
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Equation (10) requires that the .signal-to-echo ratio be large compared 
to unity. If we arbitrarily select l\/P I = 17 db, we find from (13) that 
the permissible length of waveguide is L = 18.7 miles. If we now de- 
modulate and pass the resultant signal through a baseband filter of 
bandwidth f b to obtain the FM advantage, we obtain a new baseband 
signal not unlike the original. Thus we can remodulate and proceed 
for another 18.7 miles without analytical complications, as long as the 
baseband signal-to-interchannel interference ratio remains large. It 
should be emphasized that placing a demodulate-remodulate repeater 
every 18.7 miles is done to keep out of analytical difficulties. We have 
not shown and are not implying that performance requires this crutch. 
As in the previous section, assume an FM system with f b = 10 mc and 
<r/f b = 2. We choose A/ to satisfy the inequalities (11) and (12). Pick- 
ing Af = 2 kmc yields Af/2/ = 0.018 « 1, A/ > 2 (f 6 + 4<r) = 200 mc 
and (7ttiA//2) = 25 » 1. In order to calculate interchannel interference, 
we might again multiply each echo by the factor s(nf b T) given by Fig. 
5.7 of Bennett etal. and sum the result. However, that paper gives the 
analytical expression 

s(r) = 2(a/f b )°- (tt/ & t) 4 (16) 

which is valid for the small delays of this illustration, i.e., 2irj b nT « 1 
and (27r/Vu7') 2 (a/j),)' « 1. We can integrate rather than sum. Thus in 
(10) we replace nT by t and T by dr. For the interchannel interference- 
to-signal ratio contributed by the ith trapezoid, we have 

j^ = pLLo(&or/»)y i - 2lk r (,+1)ri ^ 

jl + (1 - •-*)[*- (r/r,)]} dr 
= [ 2LLo( ' Sc ^ )VTl4r2,t ] [{6 + (i - 5)(1 - r*)J (1 + ,')■ 

- [6 + t(l -e- 2k )}i b }. 

It turns out, in this example, that the i = 0, 1 and 2 trapezoids all con- 
tribute significantly, with the i = 1 trapezoid producing the most inter- 
channel interference. In physical terms, the value of i tells us how many 
mode filter spacings the converted energy traveled in the spurious mode 
before reconverting. The ith. trapezoid represents energy which traveled 
a distance larger than i mode filters spacings but less than t + 1 spac- 
ings. In this example the i = 1 trapezoid has substantially less power 
in it than the i = trapezoid, but its larger delay more than makes up 



2776 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1963 

for this difference (see Fig. 2). The total signal-to-interchannel inter- 
ference ratio in the top channel is S/I = 66.4 db. For a total system 
length of 4000 miles, we would need 4000/18.7 = 214 repeaters and 
would obtain a signal-to-interchannel interference ratio of 66.4 - 10 
logio 214 = 43.1 db if the only contribution were due to continuous 
random mode conversion. This corresponds to 22 dba at zero level. 

This is only part of the story. The TE » spurious modes are not 
damped by helix mode niters. If it is conjectured 9 that the coupling 
coefficients, differential attenuation, etc., for these modes are the same 
for copper and helix waveguides, then the calculations and results of 
the previous section apply directly. The interchannel interference pro- 
duced by the TE02 spurious mode predominates over that from all other 
modes. Consequently, it looks as though the performance of this copper 
waveguide with mode filters is comparable with that of the helix wave- 
guide. The imposed requirement of demodulate-remodulate repeaters, 
however, is certainly a disadvantage. One must determine whether or 
not they are really necessary before any definite comparison can be 
made. 

III. MATHEMATICAL DEVELOPMENT 

Our approach is to represent the effects of continuous random mode 
conversion by an infinite number of uncorrelated echoes. These echoes 
follow directly from the expansion of the reconversion transfer function 
in a Fourier series over the frequency interval of interest. The coupled 
line equations serve as our starting point. They are 

dl (z)/dz = -l\h(z) + jc(z)I x (z) 

dl x (z)/dz = +jc(z)I (z) - Yj x {z) 

where 

7o(z), /»(«) = complex wave amplitude of the TE i mode and spuri- 
ous mode, respectively, 

r , T z = deterministic propagation constants of the TE i and spuri- 
ous modes, respectively, 

c(z) = random coupling coefficient, assumed Gaussian and pure real 
or pure imaginary, and 

z = distance along the waveguide. 

If one imposes the initial conditions 7 (0) = 1 and /*(0) = 0, then 
I (z) is the transfer function of the waveguide. This may be factored 
into two components. One component, e~ T ° l , is deterministic and repre- 
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sents the attenuation and phase which would be present in the absence 
of random mode conversion. Since in theory it is completely predictable 
and equalizable, we shall not consider this component further. The other 
component is a random variable and represents the attenuation and 
phase which is contributed by random mode conversion. This is by defi- 
nition the reconversion transfer function and is given by 

G Q (z) = e +r °%(z). (19) 

Rowe 7 shows that it may be approximated directly by use of Heard's 
method of successive approximations, or with improved accuracy it 
may be obtained from the reconversion propagation function A(z) 
defined by the relation G (z) = e _A(2) , where A(z) is approximated by 
Picard's method. We use Rowe's second approximation throughout this 
paper 

a(z) = f c ( s ) £ +Ars d S f c(t)r* Tt dt. 

Jo Jo 

The reconversion propagation function A (2) may be written as 
A(z, £, &A), where £ = &p/2ir = (ft, - 0*)/2r, AA = A - A x and 
To = A + jfto , T s = A x + JPx ■ Notice that A.4 will always be negative 
since A x is positive and larger than A Q . A is nonstationary when con- 
sidered as a function of z. Under certain conditions, however, A is ap- 
proximately stationary, when considered as a function of £, over the 
small percentage bandwidths of interest. 2 The requirement is that the 
random coupling coefficient c(z) have a white power spectrum; that is, 
the autocorrelation function R c (u) = c(z)c(z + u) be an impulse func- 
tion or a close approximation thereof. The bar denotes an ensemble 
average. This seems to be true or at least approximately true for good 
waveguides over the moderate bandwidths required in our problem. 
Fortunately, it also is true that over the small percentage bandwidths 
in which we are interested, £ is approximately proportional to frequency;* 
that is, 

£ s £0 + *' ^ m (/o - /') (20) 

where £ and / are measured at the band center, and £ = £ + £' and 
S — jo ~\- j ■ m = {v x — vo)/-c where v and v x are the mode cutoff 
factors of the TE i and spurious modes, respectively, evaluated at the 
band center frequency f , and c is the velocity of light in free space. 

* From Ref. 2, Appendix A, £ = (//2c) [(1 - y 2 )* - (1 - v/)*] = (//4c) (v** - 
t' 2 ) = constant// if v 2 « 1 and vj « 1. Letting / = /„ + /' and | = £ + £' we 
have £' + £„ = (constant//,,) (1 - /'// ) = m (f - /'). 
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Thus we may determine the mean-square value of the coefficients of 
the Fourier series expansion of the reconversion propagation function A, 
when A is considered as a function of frequency, as a limiting case of 
the power spectrum of A, when A is considered as a function of £. The 
functional form of A depends upon whether £ or / is the variable. For 
this reason we shall use a subscript 1 when £ is the variable, and a sub- 
script 2 when / is the variable. 

Formulas for the power spectrum Si(u) of the random variable 
Ai(£) are developed for the case oi AA = by Rowe and Warters, 2 
for AA = constant by Young, 8 and for the case of periodically placed 
mode niters later in this paper. 

3.1 Transfer Function and Equivalent Echoes 

The complex reconversion propagation function A 2 (/) may be ex- 
panded in a Fourier series over the interval A/ centered at f 

M.(f) = n X>ne*'" /r (2D 

where A*(/) = A,(/) for /„ - A//2 ^ / ^ /„ + A//2. T = 1/A/ and 
the subscript s reminds us that we have a periodic series approximation. 
If all of the significant energy in the FM wave, or any other wave for 
that matter, lies within A/, then the periodic series (21) may be used 
for calculations. The n ^ terms represent fluctuation components. 
If their sum is always small compared to unity with large probability, 
then the reconversion transfer function G (f) = exp [ — A 2 (/)] can be 
approximated by the series 

0,(f) = exp [- ^E «» '»""] * <""" [l - Jl «. «""" T ] ( 22 ) 

where the prime on the summation indicates that the n = term is 
omitted. The impulse response of this periodic transfer function, i.e., 
its inverse Fourier transform, is 

MO = *~"° |"*(0 - S' a n 8(t + wT)J . (23) 

Thus we have an imdistorted signal plus a series of small echoes. 
The coefficients a n are given by 

r / +A//2 

a„ = (1/A/) A 2 (/)r i2 *" //A/ df. (24) 

J/ -A//2 
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Replacing / by /' + f we have 

/•+A//2 

a n = (l/Af)r*"" oM/ / A 2 (/' +/ )6- j2Tn/ /A/ d/'. (25) 

•'-A//2 



Now, over the small percentage bandwidths in which we are interested, 
the differential phase constant is proportional to frequency: i.e., £ = £' + 
£o — ^(/o — /'). Since A2(/) = Ai(£), we may rewrite (25) as 



[mAf J-mA/12 



_ -j2Tnf IHf 

\ — —. i jiivc t owe u,c 

(26) 

where 6 n is the coefficient of the nth term of the Fourier series expansion 
of Ai(£' + £ ) over the interval -A£/2 ^ £' ^ +A£/2. For Af much 
larger than the correlation interval of Ai(|), we may approximate the 
expected value of | b n 2 1, and hence | a 2 1, in terms of the power spectrum 
of Ai(£). The result is 



\an*\ - (r/lrnDS^C-nT/w). (27) 

The imaginary part of the complex random variable Ai(£) is the 
Hilbert transform 2,8 of the real part, so the power spectrum Sx^u) 
vanishes for negative u. Thus, the coefficients a n are nonzero only when 
the sign of n opposes that of m. Negative n indicates that the echoes lag 
the signal and vice versa. Since the interchannel interference is inde- 
pendent of whether the echoes lead or lag, we may drop the minus sign 
in (27), use the absolute value of m, and restrict n to positive values. 

3.2 Uniform Differential Attenuation 

The power spectrum of the real part of the reconversion propagation 
constant Ai(£) is given by Rowe and Warters, for the case of zero dif- 
ferential attenuation and a Gaussian random coupling coefficient, as 

\(i)[8(u)L*S* + S C *[L - | u |]], \u\£L 
S(u) = , (28) 

[0, otherwise. 

Since the real and imaginary parts of Ai(£) are Hilbert transforms, the 
power spectrum of Ai(£) itself is given by* 

i(i)L 2 S c 2 8(u) + S 2 [L - u], (0-) ^ u ^ L 
SlM = (29) 

[0, otherwise. 

* Discussions of (lie autocorrelation function and power spectrum of complex 
random variables, whose real and imaginary parts are Hilbert transforms, are 
given by Dugundji, 10 Zakai, 11 and Deutsch. 12 
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Young 8 demonstrates that if the differential attenuation is nonzero 
but is a constant equal to AA, the real and imaginary parts of Aj.(£) 
are again Hilbert transforms and its power spectrum is found by multi- 
plying that for AA = by 6 - ||AAM . If | AA | L » 1 we have 

[(tfLV4)l(«) + S e *Le-* iAA]a , (0-) ^ u < oo 

&,(«) - • , . (30) 

0, otherwise. 

Equation (30) can also be obtained as a limiting case of the results of 

the next section. 

It is often convenient to work in terms of the reconversion echo-to- 
signal ratio P x /P , defined as the ratio of the power in the echoes in 
(23) to that in the signal. This is given by 

/Vft = E riFl = £ (T/\m\)S Al (nT/\m\) 

71-1 n=l 

-» [ M (l/\m\)S Al (r/\?n\)dr 
= S?L/2 \AA\. 

3.3 Mode Fillers 

Consider a long section of waveguide of length L. Assume ideal mode 
filters are placed at L , 2L , • • • , NL = L. Let each mode filter be of 
zero length and produce k nepers of differential attenuation and a 
constant differential phase shift of 6 radians, with no other side effects. 
If the differential attenuation of the length of waveguide L between 
each mode filter is negligible in comparison with that produced by a 
mode filter, we may approximate the differential propagation function 
of the waveguide by the relation 

N 

Ar(z) = j2w£ + (-/c + j6) £ 8(z - nU) (32) 

where 2x£ = A/3 = differential phase constant of the waveguide alone. 
In the coupled line equations (18) the propagation factors T and 
T x are constants. In the present case, however, they are functions of 
position along the waveguide, i.e., functions of z. Rowe has shown in 
unpublished work that in (18) we may replace 



Toz by 70(2) = [ r (5) dS 
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and 

T x z by y x (z) = f T X (S) dS, 

and use Picard's method of successive approximations to obtain the 
approximate solution 

A(L) = [ c(t)e AyU) dt [ c(S)f Ay{a) dS (33) 

"0 •'0 

The derivation and result are exactly the same as Rowe's published 
derivation 7 for constant T's with the exception of the substitution men- 
tioned. Substituting (32) into (33) we have 

A(NL ) 

c n ' jo r a i 

= J c(l) exp j2v& - (k - jO) E u(t - nL ) dt (34) 

. jf c(S) exp [-j2t£S + (k - jB) g u(S - IL )~\ dS , 

where u(x) = < . 

(0 otherwise. 

Simplifying and rearranging, we have 

AT-l n-1 r (n+l)L 

A(iVLo) = ZZ .<-* ***-" I c (/y' 2 "' dt 

„=0 /= •'n'-O 

•/ cWr^dS (35) 

J //.o 

AT-X r (n+l)/,„ /•' 

+ E / cO)*** 11 d< / c(S)e- i2 " S dS 

which after more rearranging and simplification reduces to 
A(ATLo) 

= E E ,<—«><»-*-***•» f t0 c (l + nL n )e^' r// 

n=0 i=0 •'0 

• f'° c(S + lU)f*** a dS (36) 

•'n 

+ E f ° c(/ + nLoV*** 1 d/ /" c(S + nU)e~ flTiS dS. 

n=n ^n Jn 
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If the random coupling coefficient c(z) has a white power spectrum, 
the reconversion propagation function A is an approximately stationary 
random variable 2 when considered as a function of the differential phase 
variable £. We shall now calculate the autocorrelation function Rk t (^) 
of A. Since c{z) has a white power spectrum of density S c , it has an 
autocorrelation function which is an impulse, i.e. 

Rett) = S c 8(£). (37) 

Thus, each term in the series (36) is uncorrected with all others. The 
autocorrelation function of the sum is the sum of the autocorrelation 
functions of each individual term. 

The autocorrelation function of a typical term 

" c(l + nU)^ 1 dt (38) 



f 



is 



R({) = StZ , (?!^). (39) 

In view of the fact that I < n in the double series of (36), the two 
integrals in each term are uncorrelated. The autocorrelation function 
of the double series is then 



AT-l n-l 



^ ^ e 2«-„)(*->^0) R *(fi = £ {N _ n)e -2n(*-^L ) R 2 { ^ (4Q) 

since we have N - \,n - I = I terms; N - 2, n - I = 2 terms; • • • ; 
l } n — l = N— I terms. 

The autocorrelation function of each term in the single series of (36) 
is the same as the autocorrelation function for a single piece of waveguide 
of length L Q . Thus, for c(z) Gaussian, it is the inverse Fourier transform 
of the power spectrum given by (29), i.e., it is 

\R\0) +Hfi 2 U) +jR\ii)] (4D 

where R is still defined by (39) and the circumflex stands for the Hilbert 
transform. Keeping in mind the fact that the dc terms, i.e., R(0) terms, 
add in phase, we conclude 

RaM) = (N 2 /4)R\0) + (N/2)[R\tt) +0(1)] 

a--i (42) 

+ R\Z) Z (N - n)e- 2nlk - MLo) . 

n=l 
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In order to find the reconversion echo-to-signal ratio, we could inte- 
grate over the power spectrum of Ai(£) as in (31). However, the al- 
ternative method of letting £ = in the autocorrelation function is 
easier. Note that the 0+ limit on the integral in (31) corresponds to 
dropping the constant term, i.e., the (j)P 2 (0) term in (42). Thus 
the reconversion echo-to-signal ratio becomes 

Px/Po = (N/2)R\0) + £ (N - n)r 2nk R\0). 

n=l 

If N » 1 and t~ 2Nk « 1, then 

Px/Po = ^f^ 2 coth k (43) 

where L = JVX . 

The power spectrum of Ai(£) is the inverse Fourier transform of 
(42) and is given by 

8k x (u) = (i)N 2 Lo 2 S c 2 8(u) 

N-i (44) 

+ 8* £ (N - <)« ((!■ - I « - iLo I)) 

1=0 

for u ^ (0—), where 

[x if x > 

[0 if a: ^ 0. 

If N » 1 and e~ 2Nk « 1, then the terms diminish quite rapidly. Thus, 
for i « N we have the simpler expression 

(45) 
+ NL S 2 e~ 2ik [1 + (1 - e~ 2k )(i - (u/L ))} 

where (—0) ^ i ^ u/L ^ i + 1. A rough sketch of this function is 
shown in Fig. 3. The sketch consists of an impulse at the origin plus a 
straight-line approximation to an exponential. 

It is often convenient to have an expression for the ratio of the re- 
conversion echoes under the ith straight line, i.e., from u = iL to 
u = (i + 1)L Q , to the total. This is given by 

P X JP Z = (P Il ./P )/(P,/Po) = (1 - €~ 2 V 2, \ (46) 

If we allow the differential attenuation k per mode filter and the spac- 
ing L between mode filters both to decrease without limit such that 
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their ratio remains finite, (43) reduces to (31), and (45) reduces to 
(30), where | A A | = k/L a . 

3.4 Interchannel Interference 
Let the angle-modulated wave 

ei (t) m «*">'+*<« (47) 

be applied to a filter having an impulse response given by (23), i.e., by 
MO = e~ ao W) ~ & aJU + nT)]. (48) 

<p(t) is assumed to represent stationary random Gaussian noise. The 
output wave is found by convolution and is 

eo (t) = e -*[ e **«+*« _ £' ane /»o«+nr)+,W+„r)j_ (49) 

n — -oo 

From (26) we have a n = e~ ina " T b ±n = e ien ~ jna ° T r a , where r„ and 0„ 
are real. If there are a large number of sample points 2 contained in the 
small percentage bandwidths of interest, then to a very good approxi- 
mation the n 's will be uniformly distributed in the interval (0, 27r) 
and the 6»'s will be uncorrected. Substituting, (49) reduces to 

e Q (t) = f%^ot+i*w __ g'r w l *'i*+*»+*< l +-«]. (50) 

W BO 

The output consists of the original signal plus an infinite number of 
small, uncorrelated echoes. If the total power contained in the echoes 
is much smaller than that contained in the signal, we may follow a 
procedure identical to that of Bennett et al. 5 to determine the baseband 
interchannel interference in the top channel as a function of time which 
results from each of the above echoes. We then determine the autocor- 
relation function of the interference and average over r„ and 0„ . Since 
the 0„'s are uniformly distributed, the interference caused by each 
echo adds on a power or mean-square basis. In their work, Bennett 
et al. carried along a deterministic term pT which would be replaced by 
n in our case. This term appears in the final result as cos 0„ , which has 
an average value of zero. The curves of Fig. 5.7 of Ref. 5 are for 
cos pT = and so apply to our case directly. 
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IV. CONCLUSIONS 



We have presented an analysis of the effects of continuous random 
mode conversion on an angle modulated wave. The modulating signal 
was assumed to be a band of random Gaussian noise which is used to 
simulate a single sideband frequency division multiplex. A technique 
was presented for determining the expected value of the interchannel 
interference noise which appears in the top channel at baseband. Ex- 
amples using typical state-of-the-art data demonstrated that FM ap- 
pears to be an attractive modulation scheme. Copper waveguide, with 
and without mode filters, and helix waveguide were considered. 

It should be pointed out that the techniques and the results presented 
here should be applicable to other multimode and multipath situations 
virtually unchanged. The Fourier series characterization of the trans- 
mission channel should also be of use in the study of other modulation 
schemes. 

Our work contains two key, but reasonable, assumptions which one 
must keep in mind: 

(i) The perturbation technique upon which our work is based re- 
quires that the reconverted energy consist primarily of energy which 
has been converted and reconverted once and only once. 27 

(it) The random coupling coefficient is assumed to be Gaussian and to 
have a power spectrum which is white, or uniform, over a bandwidth 
which is large compared to the bandwidth "occupied" by our signal. 
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